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ry^ , ABSTRACT 

Oh' Context. Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs 

transport magnetic flux and helicity away from the Sun. These structures are not stationary but 

5— i 1 

feature temporal evolution as they propagate in the solar wind. Simplified analytical models are 

CZ> . 

frequently used for the description of MCs, and fit certain observational data well. 
Aims. The goal of the present study is to investigate numerically the validity of an analyt- 
ical model which is widely used for the description of MCs, and to determine under which 
f*^ ' conditions this model's implied assumptions cease to be valid. 

Methods. A numerical approach is applied. Analytical solutions that have been derived in previ- 
£f} • ous studies are implemented in a 3-D magnetohydrodynamic simulation code as initial condi- 

. tions. 

Results. Initially, the analytical model represents the main observational features of the MCs. 
, However, these characteristics prevail only if the structure moves with a velocity close to the ve- 

locity of the background flow. In this case an MCs evolution can quite accurately be described 
using an analytic, self-similar approach. The dynamics of the magnetic structures which move 
5^ ' with a velocity significantly above or below that of the velocity of the solar wind is investigated 

in detail. Besides the standard case in which MCs only expand and propagate in the solar wind, 
the case of an MC rotating around its axis of symmetry is also considered, and the resulting in- 
fluence on the MCs dynamics is studied. 

Conclusions. A comparison of the numerical results with observational data indicates reasonable 
agreement especially for the intermediate case, in which the MCs bulk velocity and the veloc- 
ity of the background flow are equal. In this particular case, analytical solutions obtained on the 
basis of a self-similar approach indeed describe the MCs evolution quite accurately. In general, 
however, numerical simulations are necessary to investigate the evolution as a function of a wide 
range of the parameters which define the initial conditions. 

Key words. Sun: Coronal Mass Ejections, Sun: solar wind, Magnetohydrodynamics (MHD) 



2 Dalakishvili et al.: Magnetic clouds in the solar wind 

1. Introduction 



It is well-known that coronal mass ejections (CMEs) are one of the most significant forms of so- 
lar activity. They carry enormous masses of plasma, threaded by a magnetic field, away into the 
interplanetary medium. Further away from the Sun, these large-scale, dynamical plasma structures 
are commonly called interpla netary coronal mass ejections (ICMEs). A magnetic cloud (MC) is a 



speci fic type of ICME (see e.g. lBurlagall995MWimmer-Schweingruber et al. 2006 



magn 
alb 



Demoulin et al 



2008) and can be considered as a magnetically isolated structure moving in the solar wind. Three 
features of such magnetic structures - an enhanced magnetic field, the rotation of this magnetic fi eld 
and the low proton temperatures - are selected as bona fide signatures of MCs ( Buriagalll995 ). In 
situ observations of these physical properties of MCs are considered as important prerequisites 
towards a prediction of the geophysical effectiveness of their interaction with the Earth's magneto- 
sphere, i.e. for space weather forecasts and related issues. 

Different models for the structures of MCs have been proposed. While there is no general agree- 
ment about their large-scale structure, the local structure of MCs is commonly considered in the 



form of cylindrically symmetric force-free configurations dBurlag a 1988 



1991 



Demoulin & Dasso 



2009). It is often suggested that the ends of MCs connect to the solar surface while, according t o 



some other models, MCs are described as tori dLacoste & Ouweh and 2006; 



Romashets et al. 



2007). 



These models can be useful for capturing particular features of MCs, such as the curvature of an 
MCs axis. In a number of studies, MCs are considered as force-free, static, axially s ymmetric rigid 
flux ropes and their magnetic field is con structed on the basis of Lundquist's model ( Burlagalll988 : 
Lepping et all 1990 : Farrugia et al.lll993 ). None of these analytical studies do consider interactions 
of MCs with the ambient environment. Observations show, however, that MCs do not stay static 
but expand while propagating in the solar wind, they keep ex panding well beyond 1 AU, and they 



dynamically interact with the background solar wind flow (IBurlaga 



1998 



Demoulin et al. 2008 



Demoulin & D assollW )). 



1991 



Bothmer & Schwenn 



Analytical models describing the features of MCs are used with remarkable success by a num- 
ber of experts in order to describe certain observational data. Derivations of the analytical expres- 



a self-similarly evolving cylindrical structure ( 


Dalakishvili et al. 


2011; 


Demoulin & Dasso 


2009; 


Nakwacki et al. 


2008 


Farrugia et al. 


1995a. 


1993[ 


Leaping et al. 


1990 


BurlagJl988l). We find it 



appropriate and necessary to numerically investigate the validity of such analytical models. For 
this purpose, it is seems reasonable to start out with a comparably simple model based on ana- 
lytical considerations, and then to successively relax these simplifications towards more realistic 
configurations. This procedure is then able to reveal the extent to which simplifications such as the 
assumption of self-similarity maintain their validity in more complex settings. 

Se veral numerical studies (Ide Sterck & Poedtsll999 



2000; 



Manchester et al 



2004 



Chane et al 



2006; 



Dalakishvili et al 



2009) have been performed to explore character of magnetized plasma 



flows near MCs by t reating them as sup erconductors, i.e. the magnetic field does not penetrate the 
cylindrical structure. 



Owens et al. 



( 2005) have revealed that the region in front of an interplanetary 
MC has a rather complicated structure. In order to gain more insight into the characteristics of this 
region and the physical phenomena within it, and in order to better understand the geo-effectiveness 
of MCs, a refined investigation of MC dynamics is required. 
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For instance, it is t o be expected (and has in fact been demonstrated numerically by, e.g. 



Odstrcil & Pizzol d 19991) that variations in the structure of the background wind can severely dis- 



tort the initially simple geometry of a MC. Furthermore. iRilev & Crookerl (120041) have used MHD 
simulations as well as kinematic arguments to show that MCs tend to flatten out as they propagate 
outwards, but still stress the paramount significance of force-free field models for the interpretation 
of MC observations. It is thus of vital importance to establish under which conditions (if any) these 
simple force-free models continue to be applicable. The need to eliminate all secondary effects to 
access the true net effect of a perturbing magnetic cloud on the ambient solar wind justifies our 
deliberate choice of both a cylindrical MC geometry and an unstructured background flow (see 
Section 3) as a first step towards this goal. 

The exist ence of analytical model s properly describing the evolution of MCs is very valuable 
for the field. 



Dalakishvili et al 



(2011) numerically investigated the evolution of self-similar ana- 
lytical solutions and showed that the solutions maintain a self-similar structure for a comparatively 
long time. However, in this 1-D study the entire structure was considered to be cylindrically sym- 
metric, and it was assumed that the background solar wind could be described on the basis of a 
self-similar approach. In the present study we employ a 3-D code and implement a more general 
background flow. 



and the dynamics of CMEs ( 


V 


an der Hoist et al. 


2005 


2007: 


Jacobs et al. 


2006: 


Kleimann et al. 


2005 


2009 


Aschwanden et al. 


200alRilev et al. 


2008 


Amari et al. 


2010). 


These authors studied 



the initiation and propagation of solar eruptions in the heliosphere, and followed the evolution of 
magnetic structures in a 3-D setting. The numerical solutions show that far from the Sun, the steady- 
state solar wind can be characterized by a radial velocity and, to some approximation, by a radial 
magnetic field. Therefore, in our simulations we employed a radial flow with a radial magnetic field 
as an initial background in order to study the evolution of MCs at larger heliocentric distances. The 
diameter of the MC is only 0.2 AU at 1 AU, so we could assume the physical characteristics of the 
ambient solar wind not to change significantly in the region where it interacts with the MC. 

Additionally, since a rotation of MCs about its axis of sym metry is deduced from observations 



(Burlaga 



1995 



Farrugia et al. 



1995b; 



Klein & Burlagal ll982l) . we found it necessary to start an 



investigation of the dynamics of MCs taking into account this rotation. This is possible by ana- 
lytically formulating an appropriate initial set-up. As stated above, many groups have simulated 
MCs, and various tests and comparis ons of numer i cal models w ith each other and with measure- 
ments have been performed, see, e.g.[ 



Vandas et al 



(2010, 



2009) and reference therein. It appears, 



however, that comprehensive, systematic comparisons of analytical models of MCs with their 
full three-dimensional numeric al simulations have not be en made - there are only a few studies 
into this direction. For e xample, 



Vandas & Odstrcil 



(2000) have compared the analytical model by 
( 1995) with two-dimensional numerical simulations and found very good agree- 
ment. Although th e analy tical model served mainly as a test case for the numerical simulations, 



Osherovich et al 



Vandas & Odstrcil (2000) came to the conclusion that the analytically determined asymptotic be- 



havior of the magnetic field on the axis of an expanding flux tube has a rather general validity. 

Oth er authors have co mpared numerical and analytical results of certain aspects of MC physics, 
see e.g. lXiong et alj d2006l) . who studied the geo-effectiveness of so-call ed shock-overtaking MCs 
with a 2. 5 -dimensional model. More recently. iTaubenschuss et alJ (120101) have also performed 2.5- 
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dimensional simulati ons and compared th em regarding the expansion speed to some extent with 



analytical findings bv lOwens et al.l (120051) . resulting in average agreement. Interestingly, these au- 
thors also claim that their simulations reveal that the force-free configuration for MCs seems to be 
conserved very well, at least when averaging over the entire cross section. 

With the present analysis, we compare a recently developed analytical model with three- 
dimensional numerical simulations. The remainder of this paper is organized as follows: in 
Section 2 we describe the analytical and in Section 3 the numerical model, including the imple- 
mented initial and boundary conditions. In Section 4 we describe details of the simulation setup. In 
Section 5, the results of our simulations are presented and discussed, and in Section 6 we summa- 
rize the results of our work and discuss the conclusions as well as indicate future plans. 



2. Analytical background 



Dalakishvili et al. 



(2011) 



In this section we briefly summarize the analytical model developed in 
To perform an analytic study of the dynamics of magnetic clouds, we have to start from the full set 
of ideal MHD equations 



V-B = 0, (1) 

3,B = Vx[Fx B] (2) 

d t Q + V • (gV) = (3) 

g[d t + (V-V)]V = fx~\VxB)xB-Vp (4) 

In these equations, p denotes the thermal plasma pressure, g is the mass density, V is the veloc- 
ity field, B denotes the magnetic field, and p is the permeability of free space. 
The problem is considered in the frame of the MC and in cylindrical coordinates centered on the 
MC, i.e. with a longitudinal axis z that coincides with the MCs axis. In a number of previous stud- 
ies, the MCs were considered as cylindrical magnetic structures, characterized by axial symmetry. 
In the present consideration, both symmetry along the z axis (d z = 0) and the azimuthal direction 
(8^ = 0) are assumed. The axially symmetric magnetic field can then be expressed as 



B ee [ 0, Bp, B z ] (5) 

where B v = B v (r, t) and B z — B z (r, t). We note that this representation satisfies the solenoidal 
condition ([TJ. 

The self-similar approach, adopted here, implies that the temporal evolution of the physical func- 
tions is controlled by the self-similarity variable 



(D(f) 



(6) 



where <D(f) is a function of time. In analogy to 
equations in the following form: 



Low! (11982), let us search solutions of the MHD 
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B v = ^Q^) (7a) 

B z = O ff e z © (7b) 

q = <D«p© (7c) 

P = ^PiZ) (7d) 

One can see that the type of solutions introduced by Eqs. ( !7aH7db evolve self-similarly and 
are characterized by a particular time-scaling. Here Q v , Q z , p, and p are functions of the self- 
similar variable and <!>*, <D°", <&", and show the time scaling of the azimuthal and longitudinal 
components of the magnetic field, the plasma density, and the plasma pressure, respectively. 
We consider both a radial and a longitudinal expansion of the MC, but no motion in the azimuthal 
direction. In this case, the Eulerian velocity field of the plasma, V, can be expressed as 



V = [V r , 0, v z ] , 



(8) 



where we assume that the radial component of the velocity V r = V r (r, t), and the z component 
V z - V z (z, t), i.e. we assume that the MC maintains its cylin drical shape during its ev olution. 
The solutions readily follow from the derived equations (see Dalakishvili et akkoi 1 ), yielding 



V = 

e = 
p 



1 + kt 

Q 

(1 + ktf 
P 



[r,0,z] 



(9a) 
(9b) 
(9c) 



(1 +kt) 4 

where g and p are arbitrary functions of £ = rj <b. The components of the magnetic field read 



B 



Bo 



(1 + kt) 2 



0, /, 



i — 



M- 



(10) 



{ fo(i + kt) j \r (l + kt) ) 
where Jq(x) and J\{x) are the Bessel functions of the first kind while Bq and fo are constants. 
Here £ is a constant parameter and 1/k has unit of time and characterizes the MC's rate of radial 
expansion. 



3. Description of the model 

3.1. Coordinate systems 

It is assumed that an MC is initially a cylindrical structure placed in the radial solar wind flow. The 
MC's initial bulk velocity is perpendicular to its axis of symmetry. Hereafter, in order to formulate 
the initial conditions with particular expressions in a compact way, we introduce local and global 
coordinate systems. The global coordinate system which we use is a spherical one {R, §, tp) centered 
on the Sun. The polar axis of this system coincides with the solar magnetic axis (the z axis), and 
the azimuthal angle tp is counted from the x axis, which is directed from the Sun to the Earth. 
The local, cylindrical coordinate system (r, if, z) is related to a cylindrical magnetic cloud: the 2 
axis coincides with the axis of the MC and is perpendicular to the (x, z) plane, such that (x, z, z) 
defines a right-handed Cartesian coordinate system. The azimuthal angle of the local system <p is 
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counted from the z axis of the global system. Fig. QT| shows a sketch of the computational volume 
and the respective coordinate axes. After fixing the location of the cylinder in the global system, 
we can define functional relations between the coordinates of these two systems and transform 
vector components from one system to another. Later, the notation C r $~ z will indicate components 
of a vector C in the local coordinate system, while Cr,# ^ will denote the components of the same 
vector in the global coordinate system. 



3.2. Initial and boundary conditions 

The physical quantities are normalized as follows. The unit length Lq — 7 ■ 10 s km is equal to 
the Solar radius. The unit magnetic field Bo = 3 nT, unit number density no = 10 cm -3 , and mass 
density po = m p no = 1.1 ■ 10~ 20 kg/m 3 (where m p is the proton mass) approximate their respective 
values at 1 AU. The speed is then normalized to Voa = Bo/ y^ pp _ 20.5 km/s (i .e. of the order of 

d2009l) ) and, finally, the 



Vandasetal 



the MC's expansion velocity in the local frame, according to 
unit of time is then given by to = Lq/Vqa - 9.5 h. As a background plasma flow at large heliocentric 
distances we consider a radial flow with a radial magnetic field. In order to ensure the stationarity 
of this background flow, we choose the following expressions for the initial background density 
and radial magnetic field: 



f?out — t?0,out I n 



flouts = fio.out -I [©(#), 0, 0]. 



(Ha) 
(lib) 



Here ©(#) is an arbitrary 

function of the polar angle. In our simulations, we considered both the case 0(i?) = 1 and 
0(i?) = cos &. The latter case describes the change of sign of the magnetic field and the existence 
of a current sheet in the equatorial plane. We found that such asymmetry in the magnetic field does 
not have a significant influence on the dynamics of the MC. Inside the MC, the force-free magnetic 
fields of Eq. ( TTOb . evaluated at t = 0, is superimposed on background radial field: 



B\n - B' m + B out 
In the local coordinate system, the former becomes 



0, /i|4 , Jo (4 



(12) 



(13) 



We note that Eq. ( fT3l ) coincides with Lundquist's solution dLundquistffl950l) . 
Here, fo is a same arbitrary constant parameter introduced in Eq. ( fTOb . In the numerical simulations, 
ro is equal to the MC's initial radius. For a determination of the constants Z?o,i n and fioout, the obser- 
yational data, which show that the plasma beta inside a MC is lo wer than outside (J3 m <K l,B ollt ~ 1) 
jBurlaga et al Jl 1 98 li iBurlagall 1 99 1 : 



3j 



Bothmer & Schwenn 1998) were taken into consideration. 



The evolution of an MC is characterized by an increase of its radial and longitudinal extent, and 
by the translational and rotational motion of the whole structure in the ambient environment. The 
velocity of matter inside the MC consists of the MC's bulk velocity (which is initially perpendicular 
to the MC's axis of symmetry), the velocity caused by an increase of its radius, and the velocity 



Dalakishvili et al.: Magnetic clouds in the solar wind 7 

of lengthening. Substituting t — in Eq. ( l9al >, we derive the expressions for the components of the 
velocity inside the MC 



(14a) 



given in the coordinate system attached to the MC. The total velocity of matter inside the MC 
is then 



^ MC + K 



(14b) 



The initial background velocity is prescribed to be radial and constant far away from the MC, 
and tangential to the MC's surface, i.e. the background plasma does not penetrate the MC. An 
expression for the velocity which satisfies the above-mentioned conditions is: 



V+V 



(15a) 



where 



V=(V R ,0,O) (15b) 
in the global coordinate system. It is convenient to express V in local cylindrical coordinates: 



v r = (v s + v M c,.-K)(y) 3 

V£ = (Vmc*-V*)(^) 3 
K = VMc, 2 (y) 3 ■ 



(16a) 
(16b) 
(16c) 



Here Vs - kro is the initial Lagran gian velocity of the MC 's edge. The parameter k is the con 
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Q201 lb - Initially, the density inside the MC 



stant introduced in Eq. d9al i, see also 
is uniform and half the unit density. The initial radius of the MC is 20 in normalized units, which 
corresponds to a value of 0. 1 AU, which is confirmed to be reasonable by in-situ measurements 
iBurlagalll995L 

In addition to the kinematic case described by Eqs. (114aH16c| i. we also consider the interesting case 
in which the MC rotates around its axis. In this case we change Eq. ( I14ab and Eq. d!6bb as follows: 



K# = ™ ( 17a ) 

V'z = (V M c,e-% + ur)^ . (17b) 

Here a> is the angular velocity of the MC's rigid rotation around its symmetry axis. 

We applied the following boundary conditions formulated in the global coordinate system: on 
the inner radial boundary we prescribe for the magnetic field and density the functions given by 
Eqs. ( II laHl lbb and a radially uniform velocity. On the other (opposite and other) boundaries, the 
density, the tangential components of the magnetic field, and the velocity are extrapolated, while 
for the normal components of the velocity and magnetic field we use mass and magnetic flux 
conservation conditions. One should bear in mind that the outer radial boundary conditions (whose 
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choice was motivated by simplicity) are incompatible with the cylindrical geometry prescribed 
to the MC. While in principle one could expect an influence of these boundary conditions on 
the solutions, they turn out to be negligible far from the boundary: the propagation speed of the 
magnetic structures is the local Alfven speed. Therefore, we could conclude that the boundary- 
induced disturbances will not be able to propagate inwards against the supersonic outflow. 



3.3. Model equations and their implementation in the code 



In order to study the d ynamics of MCs num erically, we use a second-order finite volume scheme 
based on the work by Kurga nov et al.l (1200 lb for the hyperbolic part of the system of equations 
(Q]|4]i, see also 



Flaig et al 



(2009) and references therein. This is a central conservative scheme 
for the solution of equations of type d,u + V • [F(u)] = 0. In our study we solve the following 
(normalized) equations: 



d,p + V • (pV) = 
d,(pV) + V • [pVV + {p + B 2 )\- BB] = 
d,B + V • (VB - BV) = 



(18a) 
(18b) 
(18c) 



The system is closed using the ideal gas 

equation of state p = pT. The employed numerical scheme require neither a Riemann solver 
nor a characteristic decomposition, a nd was extended by means of a constrained transport descrip- 



tion for the magnetic field (see, e.g., Balsara & Spicerl ll999 



Londrillo & Del Zanna 2000), which 



a2 



ensures the solenoidality of the magnetic field. This method uses the hyperbolic fluxes to compute 
the electric field components on a staggered grid. These are then used to evolve the magnetic in- 
duction, the components of which are also given on a (different) staggered grid. The stability of 
the code and its capability t o resolve s t eep gr adie nts without i ntrodu cing spurious oscillations have 
been demonstrated, e.g. by 



Kissm ann (2006|) and 



Flaig et al 



(2009) 



In the employed code, the initial magnetic field is formulated by means of a vector potential. 
The vector potential whose curl results in the magnetic field of expression Sl3i can be expressed as 



A', 



in r,ip,z 



#0,ii 



(19) 



The vector potential corresponding to a radial magnetic field can be expressed in the global 
spherical coordinates as 



t IRJKf 



. Bo out® W sin # 

0, tp, 

R Y 



(20) 



In order to ensure continuity of the vector potential across the surface of the MC, we represent 
the vector-potential of the background magnetic field as 



Aout = Ai + V/ 



(21) 



with 



/ = #0ii 



r 1 + 



(7) ) + ~ r oJi[-^)$ + h)Jo[-^\i 



A) 



'"() 



(22) 
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Table 1. Summary of conducted simulation runs, with references to the figures which display 
selected data from the respective runs. 



Run V w B 0ou , Figures 

1 15 (200 IRf Q][6] 

2 20 (200/.K) 2 0[3]|4] 

3 30 (200/.R) 2 |5][6] 

4 20 0.05 (200/i?) 2 |5] 



while the vector potential of the MC's magnetic field is given by 

A in = A' m + A out . (23) 

4. Details of the simulation setup 

In this section, we present details of the performed simulations. The computational domain 
is a segment of a sphere bounded by R e [60,300], # e [0.1,0.9]7r, and tp e [0, Q.2]n. For the 
simulations we used 120 grid cells in the radial and polar dimensions and 60 cells in the azimuthal 
direction. 

The initial velocity of the background flow is Vr = 20, corresponding to 410 km/s. The center of 
the cylindrical MC is initially located at Rq = 170 and the MC's initial radial size is ro = 20, with 
parameter ro = 20. We conducted various simulation runs, which are summarized in Table [TJ 

In runs 1 to 3, we considered cases in which the MC's bulk velocity Vo is initialized to be 
lower than, equal to, or higher than the ambient flow speed of 20. Furthermore, the (formerly 
vanishing) MC's internal rotation around its axis of symmetry was varied, as was the direction of 
the background B field (run 4). In all runs, the mass density inside the MC was p m = 0.5 and the 
density outside the MC was p out = (200 /)R 2 . In order to describe the magnetic field inside the MC, 
we set the value of Bo.j n = 5. The Lagrangian velocity of the MC's edge is Vs = kro with k = 0.05. 

In order to ease a comparison with observations as well as with analytical models, time profiles 
of density, magnetic field magnitude, and flow velocity at selected fixed positions along the MC's 
trajectory have been extracted from the simulation data and presented in Figs.[T][3]|2l and[l0] These 
"virtual observers" thus capture information which would be measured by a stationary spacecraft 
situated on the MC's trajectory as the latter sweeps over it. All three "virtual observers" have the 
same polar and azimuthal coordinates, viz. ■&' = 0.57T and <p' = 0.037T, and the respective distances 
of each of the three observers from the Sun are R\ = 196, R' 2 = 225, and R' 3 = 250. 

5. Results and discussion 

5.1. Non-rotating MCs 

In this section, we present and discuss the results of the numerical simulations described in the 
previous section. 

Figs-EE and |7] show time profiles of normalized density, velocity, total magnetic field, and 
polar component of the magnetic field for the runs 1 to 3 as given in Table [TJ as recorded by the 
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virtual observers. In Figs. [2] [4] and [8] these numerical results are compared with those obtained 
with the analytical model for the region inside the MC, i.e. the physical variables are plotted for 
the time interval from t m when the observer situated at R = 196 enters the MC until fy,-„, when 
this observer leaves the MC again. f,„ and f/,„ are calculated analytically. These comparisons reveal 
that the analytical model describes the evolution of magnetic structures best for the case in which 
the MC initially moves with the velocity of the ambient solar wind. Fig. [5]provides plots of the 
global structure of the total magnetic field, mass density, and plasma velocity in a meridional plane 
(tp = const) and on Fig.|6]are shown structures of magnetic field and mass density in the equatorial 
plane. While these figures display the results for run 2, Fig. [9] shows the global structure of the 
number density in the meridional plane for runs 1 and 3. From the results shown in Figs. [5] and 
[6] we can conclude that when the velocity of the MC is close to the velocity of the background 
flow, the MC expands smoothly, i.e. although it does not maintain an exactly cylindrical shape 
(Fig. 13 and exhibits a slightly changing axis curvature (Pig. [6j>, the plasma density inside the MC 
decreases without developing strong gradients. Fig. |9]reveals that when the MC's velocity is less 
than the velocity of the background flow, i.e. when the MC moves slower than the background solar 
wind, a denser region appears behind the MC after some time, while the MC which initially moves 
faster than the ambient solar wind is preceded by a region with a sharp (positive) density gradient. 
Observations do indeed s how that MCs moving faster than the ambient solar wind are preceded 



by density enhancements (IKlein & Burlaga 



1982). The simulation results also show that the MC 



exhibits stronger deviations from a cylindrical shape as compared to the case shown in Figs. [5] 
and[6] From a comparison of the upper left and bottom left panels of Fig. Q] (density and absolute 
value of the magnetic field), one can see that after some time the observers detect an increase of 
magnetic field and density at about the same time (see dashed and dotted lines). Such structures 
that are characterize d by sharp gradient of density can also be found in observational data, see 
(Lepni ng et al . 1997). We could conclude t hat - besides other fa ctors such as a possible overtaking 
of slow clouds by fast corotating streams dKlein & Burlagal 1982b - dynamical processes occurring 
during the interaction of an MC with the ambient solar wind could play a crucial role for the 
creation of sharp density gradients in the vicinity of an MC. We also can see that the profile of the 
magnetic field strength evolves asymmetrically, which could be explained by the fact that the MC 
expands and, while an observer crosses this structure, the magnetic field inside the MC does not 
stay constant, but de creases in time. No te that very similar features are present in real observations, 
as exemplified in ( Lepping et al.lll997 ). We can also see that due to an increase of the MC's radius, 
the front parts of the MC are observed to pr opagate with a higher velocity than its rear regions. 
This feature is also observed in real data (e.g. Nakwacki et alj|2008 ). 



From these results it is obvious that initially, in all three cases, an observer would detect the 
following main signatures of a magnetic cloud: 1) a decrease of the plasma density, 2) an increase 
of the magnetic field strength, and 3) a rotation of the magnetic field inside the MC. We see that 
during the evolution not all magnetic structures exhibit these main features. As a matter of fact, 
in the cases when the velocity of the magnetic structure differs much from the velocity of the 
background flow, the regions of lower density and higher magnetic field evolve differently in time, 
see Figs. Q]and[7] Only those structures that move with a velocity close to that of the background 
flow maintain the above-mentioned characteristic signatures. 
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5.2. Rotating MCs 

Fig. [10] shows the MC properties for the case in which the rotation of the MC around its axis of 
symmetry z, (see FigJTTJ is taken into account. In Fig. [TO] this is done for the case in which the 
background magnetic field direction coincides with the direction of the solar wind velocity and the 
MC rotates in the north-south direction (see panel a in Fig. [TZb . In Fig. [10] peaks of density and 
magnetic field one be discerned both in front of and at the rear of the MC. This can be understood as 
follows: A rotating MC generates a centrifugal force, which induces centrifugal motion of matter. 
One could thus expect an accumulation of mass at the edges of the MC and, due to the frozen-in 
condition, the magnetic field also moves and accumulates with the plasma. In the panel showing 
the polar component B,<> of the magnetic field, we can see that the rotation causes a bending of 
the ambient magnetic field lines. For instance, in the case of a north-south rotation and when 
the background magnetic field direction coincides with the direction of solar wind velocity, B$ 
becomes negative in front of the MC, (See the bottom right panel in Fig.QJJand panel a in Fig.fTZli. 



6. Summary and conclusions 



In order to study the dynamics of MCs propagating in the solar wind, we numerically implemented 
analytical expressions for the physical variables characterizing an MC as initial conditions. These 



ous studies ( 


Burlaw 


a|| 1988 


Leonine et all 1990 


Farruda et al. 1993 


1 1995a 


Nakwacki et alJbo08 


Demoulin & Dasso 


2009 
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). 


Vandas et al. 


( 20061 12OO9I) used these functions 



to interpret particular observations at a certain time. We presented results describing the evolu- 
tion of MCs for three different cases, namely, when the bulk velocity of the MC is a) less than, 
b) equal to, and c) faster than the velocity of the ambient solar wind. We found that the results of 
our numerical simulations to be in good reasonable agreement with the observations. 

It was demonstrated that the initially prescribed main signatures of MCs, namely magnetic 
fields above ambient value, a mass density lower than that in the ambient solar wind, and and a 
rotation of the magnetic field, are best maintained in the case when the MC moves with a veloc- 
ity close to the velocity of the background flow. Due to interaction with the ambient solar wind, 



initially slow CMEs are accelerated, while fast CMEs slow down (iLvnch et al 



2003). We could 



expect that after sufficient time the velocity of solar eruption will not differ much from the velocity 
of the ambient solar wind. According to observations, approximate velocities of the MCs at 1 AU 
are 400-450 km/ s ( Klein & Burlagalll982l . 



Vandas et al 



(2009) compared observational data with similar results obtained using an ana- 
lytical approach and found them to be in good agreement. In this work the authors employed the 
same functions as we used in our present study as initial conditions. Since our numerical results fit 
observations for different times, we can further conclude that in a certain case, viz. when the bulk 
velocity of the MC is close to the background flow velocity, the MC evolves nearly self-similarly 
(see also Dalakishvili et alJboi 1 ). 

We further studied cases in which an MC rotates around its axis. It was argued that the cen- 
trifugal force leads to an accumulation of matter and magnetic field at the edges of the MC. The 
rotation is also able to cause a bending of the background magnetic field lines. We see that when 
comparing the obtained results to observations, we were able to demonstrate reasonable qualitative 
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agreement. In particular did we find cases in which the analytical models advocated by a number 
of experts continue to be valid for the description of MCs during their evolution, and such cases for 
which this is not true. 

While the numerical simulations led to interesting results which compare more easily to ob- 
servations than those from analytical approaches, we have to admit that the presented model still 
contains several idealizations. First of all, we introduced a much idealized background flow, which 
was initialized by uniform radial flow velocity, a radial magnetic field, and a spherically symmetric 



field has a logical base: a number of numerical studies ( 


van der Hoist et al.l2005 


2007 


Jacobs et al. 


2006; 


Kleimann et al. 


2005 


2009; 


Aschwanden et al. 


2008: 


Rilev et al. 


2008; 


Amari et al. 


2010) 



show that when the solar wind beyond the source surface reaches a steady state, it is characterized 
by radial plasma flow and radial magnetic field without a latitudinal component, though even if 
initially a more complicated magnetic field were introduced. It might be advantageous to start the 
simulations of the solar wind from specific initial conditions and to implement a magnetic cloud 
only later, when the solar wind has reached steady state conditions. However at this stage our nu- 
merical facilities do not enable us to perform such kind of simulations. 

Second, we considered MCs as initially cylindrical and symmetric structures. It would be interest- 
ing to also study the evolution of other more complex magnetic structures, e.g. toroids, spheroids, 
ellipsoids, as well as structures which remain connected to the Sun. In our study, the MC is initially 
already located far from the Sun and it has a certain bulk velocity. It would be reasonable to study 
the self-consistent evolution of a solar eruption from the solar surface until it reaches 1 AU. 

Since the present study mainly intends to test an analytical model, features of MCs were com- 
pared to observations only for the case of non-rotating MCs. For the future, we plan to extend these 
comparative studies also to rotating MCs. Also, in forthcoming work we intend to study the dynam- 
ical evolution of MCs in different types of background flow. Furthermore, we are also interested in 
the numerical investigation of the interaction between several magnetic structures: besides the MCs 
another class of solar ejecta, namely "complex ejecta" was identified. While most of the magnetic 
clouds were associated with a single CME, complex ejecta could have had multiple sources. It was 
conjectured that som e complex ejecta were produced by the interaction of two or more halo CMEs 
JjBurlaga et alJbood) . 
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Fig. 1. Plots of normalized density, velocity, total magnetic field, and polar component of the 
magnetic field. The solid, dashed, and dotted lines show the values recorded by virtual observers 
placed at R = 196, R = 225, and R = 250, respectively. These plots correspond to run 1 (see 
Table [B, for which V = 15 and w = 0. 
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Fig. 2. Plots of density, polar component of magnetic field, and absolute value of the magnetic 
field. The solid lines correspond to the values detected by a virtual observer situated at R — 196. 
The dashed line shows values obtained using the analytical solution. The initial bulk velocity of the 
MC is Vo = 15, its initial radius is ro = 20. 




Fig. 3. The same type of data as displayed as Fig. [TJ here for run 2, when the MC's initial velocity 
is Vq = 20 (i.e. equal to that of the ambient wind). We see that the structure maintains the MC's 
signatures during its evolution fairly well. 
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4. The same plot as in Fig. |2]but here for the case in which the MC's velocity is Vo 
to that of the background wind). 



= 20 
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Fig. 5. Cuts along the meridional plane tp = 0.037T for run 2 (see Table [T), displaying snapshots 
for the normalized values of the absolute magnetic field (top row), mass density (middle row), and 
absolute velocity (bottom row) at times t — (left column), t = 2 (middle column), and t — 4 (right 
column). The MC's initial velocity is Vo = 20 (i.e. equal to that of the background wind). The MC 
expands approximately symmetrically. 
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Fig. 6. The same situation as depicted in Fig. [3] (run 2, Vo = 20, a> = 0), except a cut in the 
equatorial plane is shown. Top row corresponds to the global structure of the magnetic field and 
bottom row shows global structure of mass density. 




Fig. 7. The same type of data as displayed in Fig.Q] here for run 3, when the MC's initial velocity 
is Vo = 30 (i.e. higher than that of the background wind). 
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Fig. 8. The same plot as in Fig.[2]but here for the case in which the MC's initial velocity is Vo = 30 
and its initial radius is ro = 20. 
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Fig. 9. Density cuts along the meridional plane at time t — 2 for run 1 (left panel) and run 3 
(right panel). In the left panel the MC initially moves with velocity Vb = 15 (i.e. slower than the 
background wind). The right panel corresponds to the case in which the MC initially moves with 
velocity Vo = 30 (i.e. faster than the background wind). We see that the slow MC is followed by a 
denser region while the fast MC is preceded by a region with a sharp density gradient. 




Fig. 10. The same type of data display as Fig.Q~|forrun 4 (Vq = 20, a> = 0.05). 




Fig. 11. Sketch of the computational volume with the global, Sun-centered (x, y, z) system, and the 
MC (cylinder with z axis and exemplary field lines) in the ecliptic plane. 
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Fig. 12. Illustration of the bending of a magnetic field line due to the rotation of the flux rope. Here 
it is assumed that solar wind flows from right to left. Panels a) and d) correspond to the north-south 
rotation of an MC, while panels b) and c) demonstrate the case of south-north rotation. Due to the 
MC's rotation field lines of the background magnetic field are curved. 



